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Abstract. We report a neutron scattering study of the instantaneous spin correlations in the two- 
dimensional spin S = 5/2 square-lattice Heisenberg antiferromagnet Rb2MnF4. The measured correlation 
lengths are quantitatively described, with no adjustable parameters, by high-temperature series expansion 
results and by a theory based on the quantum self-consistent harmonic approximation. Conversely, we find 
that the data, which cover the range from about 1 to 50 lattice constants, are outside of the regime cor- 
responding to renormalized classical behavior of the quantum non-linear a model. In addition, we observe 
a crossover from Heisenberg to Ising critical behavior near the Neel temperature; this crossover is well 
described by a mean-field model with no adjustable parameters. 

PACS. 75.10.Jm quantized spin models - 75.25.+Z spin arrangements in magnetically ordered materials 
- 75.30.Gw magnetic anisotropy - 75.40.Cx static properties 



1 Introduction 

The physics of the two-dimensional square-lattice quan- 
tum Heisenberg antiferromagnet (2DSLQHA) continues 
to receive much attention. In addition to the basic in- 
terest in studying the role of quantum fluctuations in 
this high-symmetry low-dimensional system, experimen- 
tal and theoretical efforts have heightened with the dis- 
covery that the undoped parent compounds of high-T c su- 
perconductors are typically very good realizations of the 
spin 5=1/2 2DSLQHA. In particular, neutron scatter- 
ing experiments on the insulating lamellar copper oxides 
La2Cu04 and Sr2Cu02CL> have elucidated the spin 
fluctuations of the 2DSLQHA model for the extreme quan- 
tum limit of S = 1/2. The measured spin-spin correlation 
lengths of Sr2Cu02Cl2 are quantitatively described [g] by 
quantum Monte Carlo results based on the square-lattice 
nearest-neighbor Heisenberg model H|), and by a low- 
temperature theory based on the quantum non-linear a 
model developed by Chakravarty, Halperin, and Nelson 
(CHN) Q. The remarkable agreement between the ex- 
perimental data for 5 = 1/2 and the theoretical results 
of CHN, as extended by Hasenfrantz and Niedermeyer 
(HN) S, was found not to hold for higher spin values, 
S > 1/2. In particular, in the 5=1 2DSLQHA systems 
K2MF4 pi and La2Ni04 H, the measured correlation 
lengths were found to deviate significantly from the CHN- 
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HN prediction. A systematic high-temperature series ex- 
pansion study showed that the deviations increase pro- 
gressively with increasing spin values above 5 = 1/2. Re- 
cent Monte Carlo work M suggests that the explanation 
for this discrepancy may simply be that the series expan- 
sion and neutron scattering results are not in the asymp- 
totic low-temperature regime for which the CHN-HN pre- 
diction is expected to hold. 

For the 2DSLQHA, quantum spin fluctuations strongly 
renormalize the spin-stiffness and spin-wave velocity in the 
5 = 1/2 system Sr2Cu02Ci2, but there appears to be no 
fundamental change from classical behavior. For 5 = 5/2, 
the 2DSLQHA should be even less affected by quantum 
renormalization and thus should correspond more closely 
to a classical spin system. Indeed, high-temperature se- 
ries expansion studies reveal a near-agreement between 
the correlation lengths of the 5 = 5/2 2DSLQHA and 
the 5 = 5/2 ferromagnetic system, for which the ground 
state is classical Q . In addition, upon scaling the temper- 
ature by J nn S(S + 1), series expansion results indicate 
that the correlation length of the 2DSLQHA rapidly ap- 
proaches that of the classical 5 = 00 case as the spin value 
is increased progressively from 5 = 1/2 to 5 = 5/2. In or- 
der to examine further the role of quantum fluctuations 
in the models, and their dependence on the spin quantum 
number, one requires more experimental data for the cor- 
relation length of systems well described by the 2DSLQHA 
Hamiltonian with higher spin. 

We report an energy-integrating neutron scattering 
study of the 2D instantaneous spin-spin correlations in 
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the S = 5/2 material Rb 2 MnF4. This paper is organized 
as follows: Sect. ^ contains preliminary details about the 
Rb2MnF4 system and about our measurements; Sect. || 
contains our experimental results and data analysis; in 
Sect. |], we present a comparison with various theories; 
and Sect. || summarizes our results. 



2 Preliminary Details 

We study a high-quality single crystal of Rb2MnF4, which 
is, in fact, the same as that used in previous neutron stud- 
ies 10 . Rb 2 MnF 4 has the K 2 NiF 4 crystal structure, 
space group I4/mmm, with square planes of MnF 2 sep- 
arated by two intervening sheets of non-magnetic ions. 
The magnetic Mn 2+ ions (S = 5/2) form a square lat- 
tice and arc antiferromagnetically coupled to their nearest 
neighbors via super-exchange through the intervening F~ 
ions. The low-temperature in-plane lattice constants are 
a = b = 4.20 A, and the out-of-plane lattice constant is 
c = 13.77 A. Over the temperature range of the experi- 
ment, 10 K < T < 110 K, the in-plane lattice constant 
changes by only ~0.1%, so any concomitant temperature 
variation of the in-plane exchange coupling will have a neg- 
ligible effect on the spin correlations. The exchange cou- 
pling between nearest neighbor spins in adjacent planes 
is frustrated for this body-centered tetragonal spin struc- 
ture. As a result, any effective inter-planar coupling is at 
least six orders of magnitude weaker than the coupling 
within the planes, leading to the quasi- two-dimensional 
nature of the spin system. From previous neutron scatter- 
ing experiments on this and isomorphous materials, it is 
known that the critical scattering consists of purely two- 
dimensional spin fluctuations, and no spin- wave dispersion 
is observable along the c-dircction, confirming that a 2D 
spin model is appropriate |l2^ , |il||l4| . 

The physics of the Mrr^ square lattice is well de- 
scribed by the simple spin Hamiltonian 



H — J n 



<i,j> 



St ■ s. 
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where the staggered anisotropy field Hf represents the 
effect of the dipolar anisotropy which favors spin align- 
ment along the c-axis. From both NMR measurements of 
the sublattice magnetization ]l5| | and neutron scattering 
measurements of the spin- wave dispersion [|R3]jl7j] , one ob- 
tains J nn = 7.36 ± 0.10 K. The small Ising anisotropy 
in this predominately Heisenberg Hamiltonian is aj = 
211bH A I^Zj =nn JnnSj ~ 0.0047, as deduced from anti- 
ferromagnetic resonance and inelastic neutron scattering 
measurements of the low-t emp erature spin-wave gap at 
the magnetic zone center [ p^5[jlq] . A small three-dimensional 
(3D) coupling does exist as evidenced by the transition to 
3D long-range order which occurs at T N — 38.4 K. Such a 
3D coupling usually appears as a coupling between next- 
nearest-neighbor (nnn) planes. However, the weakness of 
this coupling is evident from magnetic Bragg diffraction, 
which reveals two domains with different stacking arrange- 
ments of ordered MnF 2 planes: one in which nnn planes 



are ferromagnetically aligned and another in which nnn 
planes are antiferromagnetically aligned [jl0|. This 3D or- 
der reflects primarily 2D correlations with Ising critical 
behavior, since the staggered magnetization of both stack- 
ing domains was found to follow the same power law, 
M s = (1 - T/T N ) - 16 0; the observed critical exponent 
(3 for the order parameter is much closer to the 2D Ising 
model value (3 = | than to the conventional 3D values of 
approximately ^. 

The experiments were performed on the H4M ther- 
mal neutron spectrometer at the Brookhaven High Flux 
Beam reactor. The Rb 2 MnF4 crystal was aligned with a 
(110) axis perpendicular to the scattering plane, thus hav- 
ing the magnetic zone center wave vector (5,5,0) and 
the c-axis in the scattering plane. The spectrometer was 
operated in the two-axis energy-integrating configuration 
with collimator sequence 20'-20'-sample-20'. The scatter- 
ing geometry was chosen so that outgoing neutrons were 
perpendicular to the MnF 2 planes, thus integrating over 
energy at constant in-plane momentum transfer Q 2 d [fL3| . 
In this geometry, the intensity of the detected neutrons is 
proportional to the static structure factor, 



S(Q 2 d) 



5*(Q2D, w)dw, 



(2) 



where Q 2 d is the momentum transfer within the 2D MnF 2 
sheets, and u) is the energy transfer. Here, we assume that 
any variation in the Mn 2+ form factor has a negligible 
effect in the integration over oj. Two incident neutron 
energies, Ei = 14.7 meV and Ei — 41 meV, were used 
to verify that the experiment integrates over the relevant 
dynamic fluctuations properly. Both incident energies are 
more than an order of magnitude larger than the magnetic 
energy scale, J nn , in Rb 2 MnF4. By integrating over all 
energies, we obtain information about the instantaneous 
(equal-time) spin correlations of the system. 



W 3 Experimental Results and Analysis 



We show in Figs. 1 and 2 representative two-axis scans for 
Ei = 41 meV and Ei = 14.7 meV, respectively. In order 
to extract the intrinsic peak widths and amplitudes of the 
scattering, we fit our data to the form 



S(q 2B ) = sin 2 (</>) 



S\\(0) 



?2d/ K || 



+ (1 + cos 2 (4>)) 



S±(0) 



1 + «2dA 



(3) 



convolved with the instrumental resolution function (which 
is drawn in the top panels of the figures) on top of a sloping 
background. Here, q 2 D = (Q x — h, Qy — |, 0) repre- 
sents the displacement from the center of the rod of 2D 
scattering measured in reciprocal lattice units, and <fi is 
the angle subtended by Q and the c-axis. The correlation 
length is the inverse of the Q-space width Ku±. For 
T > 1.2 Tjv, the profiles are well described by a single 2D 
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Fig. 1. Representative energy-integrating scans along the di- 
rection (H/2, H/2, L), with L chosen such that k/ || c, for 
an incident neutron energy of 41 meV. The solid lines result 
from least-squares fits of the lineshape, Eq.(3), convolved with 
the instrumental resolution. The dashed line in the top panel 
indicates the instrumental resolution. 



Lorentzian form, that is, K|| — k±. As the system is cooled 
to the immediate vicinity of XV, the scattering exhibits a 
crossover from Heisenberg to Ising behavior [[ti] . Accord- 
ingly, we find that within about 20% of XV, for T > Tm, 
one must describe the lineshape with two components: one 
2D Lorentzian with diverging £m and <S||(0), corresponding 
to the longitudinal (|| c) Ising fluctuations, and a second 
2D Lorentzian corresponding to the non-divergent trans- 
verse (_L c) spin fluctuations. For our experimental con- 
figuration, the geometrical factors sin 2 (</>) and l+cos 2 (0) 
are both close to unity; specifically, for Ei = 14.7 meV, 
sin 2 (</>) ~ 0.96 and l+cos 2 (0) ~ 1.04. The transverse scat- 
tering contribution S±(0)/ (1 + ql D / n^) is denoted by the 
dashed line in the bottom two panels of Fig. 2. 

In our analysis, we first fit the temperature depen- 
dence of the transverse scattering at temperatures be- 
low XV using a very simple model. Within the ordered 
phase, the transverse spin contribution originates from 
spin-wave scattering; therefore, k±_ is held fixed at the 
spin-wave value of n± = 0.028 reciprocal lattice units, 
while the quantity S±(0)k]_ is assumed to increase with 



Fig. 2. Representative energy-integrating scans along the di- 
rection (H/2, H/2, L), with L chosen such that k/ || c, for 
an incident neutron energy of 14.7 meV. The solid lines result 
from least-squares fits of the lineshape, Eq.(3), convolved with 
the instrumental resolution. The dashed line in the top panel 
indicates the instrumental resolution. The dashed lines in the 
bottom two panels indicate the magnetic scattering due to the 
transverse spin component. 



the Bose thermal occupation factor. The Bose factor, 
which takes into account both neutron energy gain and 
energy loss due to spin-wave excitations, is 2n(uj) + 1 
where n{w) = l/[exp(w/T) - 1] and u = 7.2 ± 0.2 K is 
the low-temperature spin-wave gap at the magnetic zone 
center |lij . Below Tjv we also observe a sharp resolution- 
limited quasi-2D Bragg component along (i, i, L) which 
grows in intensity with decreasing temperature. This 2D 
Bragg scattering probably originates from 2D sheets at 
the interfaces between the two different 3D stacking do- 
mains. The 2D Bragg scattering is assumed to have the 
same temperature dependence as the 3D order parameter. 
This simple model is found to describe the data below X/v 
very well. 

Since the focus of our experiments is on the behavior 
above TV, the main impetus for measuring and model- 
ing the scattering below TV is to allow us to estimate 
the non-critical transverse scattering in the Ising critical 
regime above TV- By definition, the transverse scattering 
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must become identical to the longitudinal scattering in 
the Heisenberg regime well above T/y. Therefore, in the 
fitting, we assume that for the transverse component, kj_ 
and S±(0)k/j_ can be simply interpolated linearly between 
the fitted values at Tjv = 38.4 K and those at 46 K, 
and that above 46 K the two Lorentzians are identical, 
that is, within the errors the system shows pure Heisen- 
berg behavior. Thus, in the final data analysis above Tjy, 
the only free fit parameters are Ky and S'ii(O), which de- 
scribe the Heisenberg behavior at high-temperatures and 
the Ising critical scattering for temperatures near TV . The 
results are shown in Figs. 3 and 4. For the fitted val- 
ues of «|| shown in Fig. 3, the plotted error bars corre- 
spond to the larger of three statistical standard deviations 
or one-tenth of the instrumental resolution. Similarly, for 
the >S||(0) data shown in Fig. 4, the plotted error bars 
are equal to three standard deviations. We note that the 
data obtained with Ei — 14.7 meV and E{ — 41 meV 
agree well with each other, thus confirming the validity 
of the quasi-elastic approximation The dash-dotted 
lines represent the interpolated temperature dependences 
of the transverse scattering parameters. We find that vary- 
ing the upper temperature limit of the interpolation of the 
transverse fluctuations around 46 K by several degrees 
does not change the results for km and Su(0) within the 
errors. Furthermore, we also obtain closely similar results 
using a very different interpolation scheme for kj_ and 
S±(0)Kj_ similar to that employed in Ref. [jT4| . The solid 
curve in Fig. 3 corresponds to a theory for the 2DSLQHA 
by Cuccoli et al. jisj| , which we plot with a modification 
incorporating the effects of the Ising anisotropy. We will 
elaborate upon this in the next Section. The dashed line 
in Fig. 3 is the prediction from high-temperature series 
expansion for the Heisenberg model with S = 5/2 fl. In 
Fig. 4, the solid line is the unmodified result from Ref. |Q , 
and the dashed line is the series expansion result for the 
S = 5/2 2DSLQHA 0. Since the neutron scattering in- 
tensity is not measured in absolute units, we choose a scale 
for the plot of S\\ (0) in Fig. 4 such that the data approach 
optimally the limit Sj(0) -> S(S + l)/3 as T -> oo. 

Previous studies |Q of the 2D antiferromagnets K2N1F4 
(S = 1) and K 2 MnF 4 (S = 5/2) show that within ~20% 
of Tjv the critical magnetic scattering follows the behavior 
predicted for the 2D Ising model, for which the exact ex- 
ponents are v—\ for the correlation length and 7=1.75 for 
the susceptibility. Our main purpose here is to study the 
2D Heisenberg regime, so the instrumental resolution was 
not optimized to investigate the narrow, rapidly diverg- 
ing peaks in the Ising critical regime. Even so, our fitted 
values for km at temperatures within 20% of T/v yield an 
exponent of v = 1.0 ± 0.1 in good agreement with the 
exact result v = 1 for the 2D Ising model. 



4 Comparison with Theory 

A low-temperature theory for the 2DSLQHA was for- 
mulated by Chakravarty, Halperin, and Nelson, in which 
they obtained the static and dynamic properties of the 
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Fig. 3. Temperature dependence of Ky , the width of the longi- 
tudinal spin component of the critical scattering. At high tem- 
peratures, km exhibits 2D Heisenberg behavior, crossing over 
to 2D Ising behavior near Tjv. The dash-dotted line represents 
the interpolated temperature dependence of ki . The solid line 
is the PQSCHA result [18] modified to include the effects of 
the Ising anisotropy, and the dashed line is the prediction from 
series expansion for the S = 5/2 2DSLQHA [7]. The plotted 
curves have no adjustable parameters 
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Fig. 4. Temperature dependence of S\\ (0), the peak intensity of 
the longitudinal critical scattering. The dash-dotted line repre- 
sents the interpolated temperature dependence of 5^(0). The 
solid line is the PQSCHA result [18] , and the dashed line is the 
prediction from series expansion for the S — 5/2 2DSLQHA 
[7]. Since the neutron scattering intensity is not measured in 
absolute units, the data are scaled to approach optimally the 
limit 5|| (0) -» S(S + l)/3 as T -> 00. 
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2DSLQHA by mapping it onto the 2D quantum non- 
linear a model Q . The 2D quantum non- linear a model is 
the simplest continuum model which reproduces the cor- 
rect spin- wave spectrum and spin- wave interactions of the 
2DSLQHA at long wavelengths and low energies f|,|ljl- 
First, CHN argued that the 2DSLQHA corresponds to the 
region of the 2D quantum non-linear a model in which 
the ground state is ordered - the renormalized classi- 
cal regime. Then, CHN used perturbative renormalization 
group arguments to derive an expression for the correla- 
tion length to two-loop order, showing a leading expo- 
nential divergence of £ versus inverse temperature. Later, 
Hasenfratz and Niedermayer employed chiral perturbation 
theory to calculate the correlation length more precisely to 
three-loop order ||. In the renormalized classical scaling 
regime, the correlation length is given by Q 



e c/a 
8 2np s 



a 2-K P s/T 



1 / T 

2 \2npt 



O 



( T 



(4) 

which we refer to as the CHN-HN formula. The parame- 
ters p s and c are the macroscopic spin-stiffness and spin- 
wave velocity of the model, respectively. For the nearest- 
neighbor 2DSLQHA, they are related to the microscopic 
parameters J n mS and the lattice constant a according 
to p s = Z p (S)S 2 J nn and c = Z c (S)2\/2aS J nn - The co- 
efficients Z p (S) and Z C (S) are quantum renormalization 
factors, which can be calculated using spin-wave theory 
(S > 1/2), series expansion (S =1/2,1), and Monte 
Carlo techniques (S = l/2)_M^§. For S = 1/2, 
the spin-wave approximation [ijj gives Z p (l/2) ~ 0.699 
and Z c (l/2) ~ 1.18, whereas for S = 5/2, the factors 
Z p (5/2) ~ 0.951 and Z c (5/2) ~ 1.03 are closer to the 
classical limit of Z p = Z c = 1. 

As mentioned above, the correlation length of the 
S = 1/2 system Sr 2 Cu0 2 Cl 2 is well described by Eq. (4) 
throughout the entire experimental temperature range. 
Subsequent Monte Carlo work § indicates that this agree- 
ment for the S = 1/2 case is, at least in part, coincidental; 
deviations exist, but they are too subtle to be discerned 
experimentally. In CHN-HN's formulation, the natural ex- 
pansion parameter for temperature is T/2irp s - For all ex- 
perimental systems, the lowest temperature at which 2D 
Heisenberg behavior can be observed is bounded by a non- 
zero temperature T/v below which there is 3D long-range 
order. For the S = 1/2 system Sr 2 Cu02Cl2, the mea- 
sured experimental temperature range is 0.16 < T /2irp s < 
0.36 0. In our present study of Rb2MnF4, which has a 
higher spin, though smaller J„„, the temperature range 
is similar, 0.14 < T/2np s < 0.40. The Monte Carlo sim- 
ulations of Beard et al. || for S = 1/2 indicate that the 
three-loop CHN-HN formula, Eq. (4), provides an accu- 
rate description of the S — 1/2 2DSLQHA for tempera- 
tures T /2tt p s 5 0.15. This barely overlaps the temperature 
range studied experimentally here. Further, for S = 5/2 
the renormalized classical scaling regime is expected for 
lower temperatures [Q. 

In Fig. 5, we plot our results for the correlation length 
versus inverse temperature. Since J nn is known from inde- 
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Fig. 5. The longitudinal correlation length £y = 1/k\\ ver- 
sus inverse temperature in reduced units. The open symbols 
are our data for Rb2MnF4. The closed circles are results for 
KFeF 4 from Ref. [21]. The dashed line is the PQSCHA [18] re- 
sults modified to include the Ising anisotropy. The dash-dotted 
line is the unmodified PQSCHA result. The solid line is the re- 
sult of series expansion for the S = 5/2 2DSLQHA [7]. The 
dotted line is the three-loop CHN-HN formula, Eq. (4) [4,5]. 
The comparisons are in absolute units, with no adjustable pa- 
rameters 



pendent experiments |15[|16| ], there are no adjustable pa- 
rameters in the comparison between theory and experi- 
ment. Also plotted in Fig. 5 are the experimental data 
from Fulton et al ^ for the 2D S = 5/2 Heisenberg 
antiferromagnet KFeF4. The KFeF4 system is less ideal 
because the nearest neighbor spins in the planes are dis- 
torted from a square configuration forming a rectangular 
lattice. However, using the average of the magnetic ex- 
change coupling along the two different in-plane lattice 
directions (J^ s — 26.9 K from neutron scattering mea- 
surements of the spin- wave dispersion |nj with zero point 
correction |l7|), we find that the agreement for the mea- 
sured correlation lengths between Rb 2 MnF4 and KFeF4 
is excellent. The reduced Ising anisotropy for KFeF4 of 
ai ~ 0.0045 is fortuitously almost identical to that for 
Rb 2 MnF 4 . The three- loop CHN-HN formula is plotted 
for S = 5/2, and it appears as a straight line on this 
logarithmic scale. We used the spin-wave theory values of 
p s — 5. 94 J and c ~ 7. 30 J a [0. It is evident that the data 
deviate significantly from the three-loop CHN-HN result 
for the quantum non-linear a model. 

Figure 5 also includes the series expansion result for 
the S = 5/2 2DSLQHA from Ref. §. There is good agree- 
ment between the series expansion prediction, with plot- 
ted correlation lengths up to 14 lattice constants, and our 
data over most of the temperature range. Systematic devi- 
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ations exist at very high temperatures, where the correla- 
tion length is small, and at low temperatures near Xjy, 
where the correlation length begins to diverge rapidly. 
At high temperatures, where the correlation length de- 
creases to one lattice constant or less, £ is not precisely 
defined. Also, at such high temperatures where the mea- 
sured peaks become broad in reciprocal space, the use of 
an Ornstein-Zernike Lorentzian lineshape may be called 
into question, and the background is more difficult to de- 
fine. Accordingly, it is not surprising that slight deviations 
exist at high temperatures. We have already pointed out 
that Ising fluctuations will dominate near T/v, so the cor- 
relation length will necessarily deviate from Heisenberg 
behavior in the immediate vicinity of TV. 

An alternative theoretical analysis of the 2DSLQHA 
has been carried out by Cuccoli et al. |Q in which they 
treat quantum fluctuations in a self-consistent Gaussian 
approximation separately from the classical contribution. 
In their approach, which they label the purely-quantum 
self-consistent harmonic approximation (PQSCHA), the 
quantum spin Hamiltonian is rewritten as an effective clas- 
sical Hamiltonian, where the temperature scale is renor- 
malized due to quantum fluctuations, and the new classi- 
cal spin length appears as 5+^. Defining the reduced tem- 
perature as t = T/{J nn (S + \) 2 }, the correlation length 
for the 2DSLQHA is then simply given by 



£(*) = £ci{t c i) with t c i = 



t 



e\ty 



(5) 



Here, £ c / is the correlation length for the corresponding 
classical 2D square-lattice Heisenberg model, and 9 4 (t) is 
a temperature renormalization parameter. The PQSCHA 
is most accurate in the limit where the quantum fluctu- 
ations are weak, and correspondingly O^jt) is near unity. 
The calculations of Cuccoli et al. 18 show that this is 
the case for 5 = 5/2 over an extended temperature range; 
further, their results show good agreement with the exist- 
ing experimental S — 1/2 and 5=1 data over the ap- 
propriate high-temperature ranges. The behavior of £ c ; is 
determined from classical Monte Carlo simulations. The 
PQSCHA result for 5 = 5/2 is plotted in Fig. 5 as the 
dash-dotted line. Similar to the series expansion result, 
there is good agreement between the PQSCHA theory and 
our data and those in KFeF4 with no adjustable param- 
eters. Again, as TV is approached, the data deviate from 
the theoretical curve because of the crossover to Ising be- 
havior. 

Keimer et al. jl| have derived a simple mean-field 
model which allows one to incorporate the effects on 
the longitudinal spin correlations of the staggered Ising 
anisotropy field as appears in Eq. (1). In this model, the 
unperturbed Heisenberg correlation length is replaced by 
the form 

g(a/,T) = — (6) 
V 1 -a/^ s ( T ) 

Using the value a/=0.0047 measured in Rb2MnF4 and 
the form for £ He is{T) given by the PQSCHA for the 
Heisenberg model, we obtain the dashed line shown in 



Fig. 5. This mean- field result for the correlation length 
nicely captures the crossover from Heisenberg to Ising 
spin correlations in both Rb2MnF4 and KFeF4. By de- 
sign, this model takes into account the growing impor- 
tance of the local anisotropy field as the spin correlation 
length grows; it does not include 2D Ising critical effects. 
Thus, the predicted divergence is mean-field-like with a 
power-law exponent of v—1/2 instead of the 2D Ising re- 
sult of v=l. This rapid mean- field divergence is also ele- 
gantly seen as the "PQSCHA-with-anisotropy" curve rises 
above our data very close to T/v. Nevertheless, this mean- 
field model with no adjustable parameters predicts TV 
to within 3% in Rb 2 MnF 4 and to within 6% in KFeF 4 . 
The accuracy of this model is highlighted when com- 
pared to the conventional mean-field theory prediction of 
T^v = 4J n „5(5 + l)/3 which is higher than the experi- 
mentally measured transition temperatures by more than 
a factor of two. 

In Fig. 6 we plot the ratio 5|| (0)/£jj versus temperature 
for Rb2MnF4. The open symbols indicate data for which 
the Ising critical fluctuations become significant, so we 
concentrate here on the closed symbols in the temperature 
regime where the behavior is predominantly Heisenberg- 
like. The low-temperature analysis of CHN-HN predicts 
that this quantity should follow the form 



5(0) 

e 



A2nM? 



T 
2irp s 



1 + C 



2np s 



O 



( T 



\ 2np 6 



(7) 

where M s is the T = staggered magnetization and A 
and C are universal constants. The same form, which is 
written above at three-loop order, can also be derived 
from a renormalization group analysis of the classical 
model |p|,p2|. In an alternative approach, Kopietz 24 ob- 
tained the leading T 2 behavior based on a Schwinger bo- 
son mean-field theory (23) . Using the values for the univer- 
sal numbers A and C in Eq.(7) obtained from Monte Carlo 
calculations on the 5 = 1/2 and 5=1 models (A — 4.5 
and C = 0.5) 25 2fJ, we plot Eq.(7) in Fig. 6. It appears 



that the low-temperature data points in the Heisenberg 
regime may overlap with the T 2 law; however, the data 
depart significantly from the predicted T 2 law at higher 
temperature. Since the overall scaling factor for the S\\ (0) 
data was determined in Fig. 4, there are no adjustable 
parameters in this plot. Also plotted are results from se- 
ries expansion § and the PQSCHA theory |Jj for the 
5 = 5/2 2DSLQHA. There is a rough overall agreement 
with both of these latter results, but systematic deviations 
clearly exist. 



5 Summary 

In summary, we find that the instantaneous spin-spin cor- 
relations in Rb2MnF4 are quantitatively described by high- 
temperature series expansion results and the PQSCHA 
theory for the 5 = 5/2 2DSLQHA with no adjustable 
parameters. There is also good absolute agreement with 
previous results for KFeF4 , which is a somewhat less ideal 
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Fig. 6. The ratio S||(0)/£y versus temperature. The dashed 
line is the PQSCHA result [18]. The dash-dotted line is the 
prediction from series expansion for the S — 5/2 2DSLQHA [7]. 
The solid line is Eq.(7) with the parameters A and C obtained 
from [26]. 



S = 5/2 2DSLQHA. Our data, which correspond to corre- 
lation lengths up to about 50 lattice constants, are not well 
described by the three-loop CHN-HN prediction based 
on the 2D quantum non-linear a model. This discrep- 
ancy almost certainly derives from the fact that the ex- 
perimental data correspond to temperatures where the 
5 = 5/2 2DSLQHA is not in the renormalized classical 
low-temperature regime described by Eq. (4), but rather 
in the classical scaling regime. Arguments by Elstner et 
al. , based on CHN's use of cutoff wavevectors for inte- 
grations over the Brillioun zone, give a crossover tempera- 
ture T cr ~ JnnS between renormalized classical behavior 
at low-temperature and classical scaling behavior at high- 
temperature. For Rb2MnF4, T cr ~ 18.4 K, which sug- 
gests that the experimental data lie in the classical scaling 
regime. In addition, the PQSCHA theory predicts the cor- 
relation length accurately in absolute units for this mate- 
rial, which is consistent with the statement that quantum 
fluctuations are not large and the system is in the classical 
scaling regime in the temperature range measured. 

As of yet, no quantum Monte Carlo results exist for 
the S = 5/2 2DSLQHA. In Monte Carlo work for S = 
1/2 PJ25| it was found that £ deviates somewhat from the 
three-loop CHN-HN result for correlation lengths less than 
~ 1000 lattice constants. However, results for the ratio 
S(0)/£, 2 f2f| with fitted universal constants show quan- 
titative agreement with the T 2 law of Eq. (7) over the 
complete temperature range. Recent Monte Carlo work 
for S = 1 [Q shows that for correlation lengths of less 
than 25 lattice constants, £ is not well described by the 
three-loop CHN-HN low-temperature formula. However, 
the S = 1 Monte Carlo data match well both the results 



from the PQSCHA theory for S = 1 and the experimen- 
tal data for La2Ni04. Also, the leading T 2 behavior for 
S(0)/£ 2 is followed by the Monte Carlo data for 5 = 1. 
This is in contrast with our results for Rb2MnF4 where we 
find a significant departure from the leading T 2 behavior 
for S(0)/£ 2 at high temperatures or, equivalently, short 
correlation lengths. 

Finally, we note that the 2D Heisenberg-Ising crossover 
behavior in both Rb2MnF4 and KFeF4 are predicted rea- 
sonably well by Eq. (6) without any adjustable parame- 
ters. However, since Eq. (6) is a mean-field result it cannot 
account for the asymptotic 2D Ising critical behavior. 

In conclusion, we now have a quite complete under- 
standing of the spin correlation length in the S — 5/2 2D 
Heisenberg antiferromagnet for the experimentally rele- 
vant temperature range. Further work is required to un- 
derstand the static structure factor peak intensity S(0) to 
the same degree. 
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AC02-98CH10886, Division of Materials Science, U.S. De- 
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